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1 . INTRODUCTION 


Buffeting flow arises when flow separation occurs on an 
airplane. The resulting flow field is highly turbulent, thus 
producing fluctuating pressures on lifting surfaces in the detached 
flow region. Boundary layer separation is perhaps the most common 
source producing buffet on most conventional configurations. 

Research in this area has been quite extensive and involved 
measurements of fluctuating pressures on models together with some 
theoretical methods to extrapolate these results to full-scale 
vehicles (see, for example. Refs. 1—3). Frequently, these pressure 
measurements are made on a conventional "rigid" model, instead of an 
aeroelastic one, because the latter may not be able to withstand 
high enough dynamic pressures to be realistic. Based on this 
consideration, several theoretical methods to use these pressure 
measurements to predict buffet response have been developed. Some 
of these methods will be reviewed later. Review of some test 
results can be found in References 4 and 5; and of theoretical 
methods, in References 6 and 7. 

Of particular interest in the present investigation is the 
buffeting caused by leading-edge vortices on slender wings. Test 
results showed that 

(1) buffeting was low before vortex breakdown and became 
severe after that (Refs. 8 and 9); 

(2) Low-frequency buffeting was more severe (Ref. 8); 
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( 3 ) 


high-frequency buffeting was caused by boundary layer 
fluctuation, and leading-edge vortices produced mainly 
low-frequency fluctuation (Ref. 10); 

(4) the results were not sensitive to Reynolds numbers (Refs. 
10 and 11), so that flight and tunnel measurements could 
be well correlated (Ref. 12); 

(5) buffeting at vortex breakdown was associated with Che wing 
response at the fundamental mode (Ref. 8). 

One conclusion from this early-day research on leading-edge 
vortices was that the buffeting induced by vortex breakdown would 
mostly be academic because a slender-wing airplane would normally 
not operate in the vortex-breakdown region of angles of attack. 
Investigation on the effect of vortex breakdown on the buffeting of 
nearby lifting surfaces, such as tails, was scarce. However, it is 
known that the vortex from the strake (or leading-edge extension, 
LEX) may reduce the buffet intensity on the wing before it bursts 
(p. 109, Ref. 7). 

In the present study, the main objective is to predict 
buffeting on vertical tails induced by LEX vortex bursting. 
Fundamental equations for structural response will first be 
derived. Existing theoretical methods for buffet prediction will be 
reviewed. The present method and some numerical results wi L 1 then 
be presented. 

In the Appendix, results of water tunnel testing of an F-18 
model are described. 





2. THEORETICAL DEVELOPMENT 


2.1 Formulation of Equations 


Structural Equations of Motion: 

Let the structural displacement, z a (x, y, t) , be expressed in 

terras of normal mode shapes, <j> n (x, y) . Then 

N 

z a (x, y, t) = l q (tH (x, y) (l) 

n=l 

where q n (t) is the so-called generalized coordinate . It can be 
shown that the structural equations of motion In forced oscillation 
in generalized coordinates can be written as (Ref. 13, pp. 131-139, 
or Ref. 14, Chapter 10): 

M n q n + Vn 2q n = Q n = ff lp E + P M J *n U ’ n)d ^ dl1 (2) 

where 

2 

M n = / / <|> mdxdy, the generalized mass 

m(x, y) = mass per unit area 

= frequency of the n cn normal mode 
Q n = the generalized force. 

The generalized force consists of two terms, one being the 
externally applied force (i.e., the p E -term) and the other being the 
force due to structural motion (i.e., the p M -term). 

The p^— term can be further decomposed in terms of the 
generalized coordinates as 


M 


N q . 

^ A Pj (x, y; u>, Mj ^ 


(3) 
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where Apj is the lifting pressure at point (x, y) on the wing caused 
by the motion of the normal mode and b Q is the reference length, 
e.g. the root seraichord. It follows that 


<Mn 


N q, N q 

I b // Ap.^d^dn = q a I fl AC HCdn 

j-1 o S w j=l o S w F j 


q N 

= b l — ) A , q , 

o b_ , nj M j 


(4) 


o j = l 

where A n j is the generalized aerodynamic force matrix and is defined 
as 


A nj = // AC p ♦a (5) 

j 

5= b o c 

(6) 

n = 

2 

In Equation (4), q^ is the dynamic pressure ( = pV^ /2). 

Equation (2) can now be written as 


M q + M u) q 
n n n n n 




N 

\ A . q . 
. , nj j 
J=1 


// p„U,n,t)<|> (5, n)d£dn 
r* n 


- V (t) 


(7) 


In the above derivation, neither structural nor viscous dampings 

have been Included • To include the former, a) - is usually replaced 
2 

with a) ( 1 + ig^) y where g n is the structural damping coefficient 
t h 

for the n mode and Is usually taken to be 0.03 if not known 

experimentally. To account for the latter, 2 CM u) q Is added to 

n n n 

the equation with £ being the damping ratio. Equation (7) becomes 
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( 8 ) 


• 9 ' 

M q + 2?M a q + M cu (l+ig)q - fcq 7 A .q. = Q E (t) 
n n n n n n n n ^°° L m v J 

j=i j j 

Structural Response to Random Excitation: 

If the excitation force Q n ^(t) is random, it may be represented 
in a Fourier integral (Chapter 14, Ref. 15), 


where 


Q n (0 - / Q n (iu>)e du> 

— 00 


Q n E (iw) - lim-L J 0 E (t)e~ iuit dt 

T*oo -x 


(9) 


( 10 ) 


The displacement q n (t) will also vary randomly, so that a Fourier 
integral representation is appropriate. 


q (t) = / q ( ioOe^^du 


( 11 ) 


Substituting Equations (9) and (11) into Equation (8), and requiring 
the relation to be valid for all t, it is obtained that 

2 2 N 
q n (i“)[- M n“ + 2i VV* + M n ^ n (1 + ig n >] - *q. I A q 

j=l J J 

= Q n E (ia») 


or. 


l { [ — M <i» + 2iM (u u + M u) 2 (l + ig ) ] 5 . - £q A . }q. 

j __ I **■ n n n n n nj ~ °o nj ^ j 


= Q n E (iai) 


n = 1 , . . . ,N 


( 12 ) 


5 



Let 


2 2 

Z . ( u>) * [ -M u> + 2iM £o) 0) + M a) ( 1 + ig ) ] 5 . 
nj n n n n n n nj 


- A . 
00 nj 


(13) 

Note that Z n j(a)) is called the complex impedance of the system; and 

its inverse, , is the so-called structural transfer function. 

To describe quantitatively a random response in a meaningful 

manner, statistical methods must be used. The most important 

quantity for this purpose is the mean square value. It is defined 

for a random function F(t) as (Ref. 15) 

X X 00 

F 2 (t) * lim / F 2 (t)dt = lim ~ / F(t) / f ( iuOe* ^durit 

X*ao —X X+°° Z 1 — T — « 


= lim ^ / f(ia>) / F( t)e la)t dtd 
X>« ~ —op —X 


to 


op op 9 

= lim ■ / f(iu>)2nf (im)du> = / lim * - du» (14) 


2T 

T-H» -CD 


— OP X 


where f is the complex conjugate of f. Define 


S('o) = lim 
X 

Equation (14) becomes 


irf(iui)f (iu») 


(15) 


F Z (t) = / S( u>)d 


OJ 


(16) 


In case the random function depends also on space coordinates, 
the definition of S(u») must be modified. For the generalized force 

til p 

of the n mode, Q n (t), it is defined as (Equation 7) 


b 



Q n E ( c ) ■ II P E (r, t) <|> n ( r)dA 


(17) 


where space coordinates (5, h) are now represented by r. The 
Fourier spectrum of Q n ^ is 


Q n E ( i(1) ) = II P E (ico) <j> n (r)dA 


(18) 


The power spectrum of the n ch generalized force is given by 


S n (o)) = lim j Q n E (iu>)Q * E (io)) 
t+°° 


llm 7 / Pp(lw) <f> (r )dA I p ( iu») <J> (r )dA 
t-m» 1 a ^ n x i n i 


T+® A 


A 

T T 

f I 

4 7r^ -T -T 


' £t {[ ^ L i " E <V WV V 


iu)(t -t ) 

e dtjdt 2 dA^dA^ 


(19) 


Let ^2 ~ t[ - x. Equation (19) can be written as 

T T 

VW r 2'“" . !. 

AA -T -T 


S n (u)) = If ♦ n (r 1 )* (r )ll.iij / / p (r t )p (r t + t) 

A A T -woo —T 1 t* ± 1 u *• 1 


— 1 OJT 

• e dt ^ dxdA^ dA ? 


OO 

^n^ r l^n^ r 2^ ^ 2 tt E l2^ r l’ r 2’ T ^ e dtdA^dA^ 


where 


1 T 

R 12 (r i* r 2 ’ T) = ^2? / T P E (r i* t l ) P E (r 2’ t L + T)dt l 

is known as the space-time correlation or cross correlation 


( 20 ) 


( 21 ) 
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function. Define 


S 12 (r 1 ’ r 2’ ~ hr f R i2 (r l’ r 2 * T)e iU>TdT (22) 

— oo 

the cross power spectral density of pressures at r^ and r 9 . It 
follows that 

S n (u>) = / ^(rj) / ♦ n < r 2 )s i2 (r i» r 2’ u,)dA 2 dA i (23) 

A A 

S ( oj) or S n (ui) is known as the power spectral density. This is 
because if F(t) were a current, the power developed by this current 
as it passed through a resistance of one ohm would be F (t). 

Returning to calculation of the total response. Equations (11) 

- (13) show that the amplitude of the motion in the n^ mode is 

00 

q n (t) = / ( [ Z( id) ] _1 {q E ( iij)) } ) n e iU>t du> (24) 

— 00 

The total displacement Is therefore 

N °® 

z (x, y, t) = l I ([Z((D )] -1 {q E ( ia») }) * (x, y)e 1U,t d(D (25) 

a , n n 

n=l -°° 

from which the Fourier spectrum of the total displacement can be 
identified as 

N -I — F 

l Clz( u») ] {Q (iui) } ) n $ n U, y) ( 26 ) 

n=l 

and the corresponding power spectrum is 

N _ N 

S (<d) = Lim - { l ([Z(«)] {Q ( iu)> }) d) (x, y ) } { £ ([Z(u>)] 1 
T-h» n=l n n n=i 


{Q E (iu) } ) n ^(x, y) } 


(27) 
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Once S w (a>) is knovm, the mean square value of displacement can 
be obtained as 


2 


z 


a 


00 



S (u))du) 
w 


(28) 


Responses in accelerations, loads, moments, and stresses, etc., 
can be similarly formulated. 

Equation (27) is difficult to simplify because of mode coupling 

through the generalized aerodynamic force matrix, A n j . If the 

aerodynamic force due to structural motion is ignored, or A » Q if 

nj 

n ^ j, then Equation (27) can be further simplified. Let 


Z nn U) " M n [_a) + 21 «“ n “ + “n^ 1 + ig n )] " 


(29) 


Equation (27) can be rewritten as 

~ E 


N Q E (ia))4 (x, y) N Q* E (iu) + (x, y) 

V M) 3 lim f { l 2 H l ** - 

Tx» n-l Z nn (a>) n=l Z ( to) 


Z nn* (U)) 

nn 


After multiplying this out, it can be obtained that 

- P _p* 0 

N Q Q $ (x, v) 

s w C«> - u- f { I ■ n ? 1 + 

T+« n =l Z Z 

nn nn 


_ p_ *p 

N N Q. Q. 

- I l 

j=1 * =l z \\ z %% 


m.i 


* r j r i- 


N <J> n (x, y) 


(30) 


,z ( B ,|2 1 l *n (t 2 >*li< r r V 

nn 1 

N N 4>j(x,y)^ 4 (x, y) 


L l 

j=l 4-1 Z..(w)Z. a (oj) A 


j*<i 


jj U 


/ ♦j c r t ) / 4*^^ r 2^ S l2 < - r l * r ? * u>)dA 9 dA 1 

(31) 
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The first series of Equation (31) represents the sura of the spectra 
of the responses in individual modes. The second series represents 
the correlation between the responses in different modes. The 
second series can be ignored if only two or three modes are present 
and their natural frequencies are widely separated (Ref. 15). 

In Reference 16, the cross power spectral density was specified 
in exponential functions with coefficients determined by experiment. 

2.2 Existing Theoretical Methods for Buffet Prediction 

All existing theoretical methods require some types of 
experimental data to work with. Sophistication of these required 
data distinguishes one method from the other. 

Cunningham and Benepe (Refs. 17 and 18): 

Pressure power spectral densities are first converted into 
pressure distributions over the wing for each frequency. The 
doublet lattice method (DLM) is then used to calculate Induced 
pressures on the tail due to downwash produced by the wing buffet 
pressures. The wing and tail pressures are used in the DLM to 
calculate the generalized aerodynamic forces. The whole equation 
(12) is used without simplification. The calculation is similar to 
that for gust response. 

B. H. K. Lee (Ref. 19): 

Again, Equation (12) is used. However, the cross correlation 
function Sp in Equation (23) is either taken to be constant over an 
aerodynamic panel or asumed to vary exponentially in space. 
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Mullans and Lemley (Ref. 20): 

The fluctuating pressure on a rigid model is again used to 
calculate the generalized aerodynamic forces. However, the 
aerodynamic forces due to wing vibration, i.e. A^-terms, are 
ignored. 

J. G. Jones (Refs. 21 and 22): 

It is assumed that each mode behaves as a single— degree— of- 
freedom system: 


Mq + 2M £u> q + M (0 ^q = Q ^(t) 
n n n n M n n n M n y n ' 


(32) 


The aerodynamic forces due to wing motion are ignored. Applying the 
Fourier transform to Equation (32), it is obtained that 


M (-U) 2 + 2i£w a) + a> 2 )q = 0 E ( i u>) 

n n n n n x n 


(33) 


Using the definition of power spectral density. Equation (15), the 


power spectral density of the response can be obtained: 
S (w) = 


s q <“> 


2 2 2 2 2 
M (-oj + 2i£u> w + w )(-u> - 2iCu) o> + u) ) 

n n n n n 


(34) 


The mean square value of q n is therefore 


q n ^ = / s (^dw 

— A ^ 


* / 


s- <«,) 

< n 


M 2 H (u»)H *(«) 
n n n 


du) 


(35) 


where 


2 2 
H (u>) = -u> + 2itu) in + <i> 

n n n 


(36) 
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The main contribution to the value of the integral in Equation (35) 
comes from the peak response at id = If S-( u>) is assumed not to 

vary appreciably in the neighborhood of o> n , it can be factored out 
of the integral in Equation (35) and the result integrated 
analytically based on the residue theorem in the theory of a complex 
variable. Results are available in Reference 23 (p. 218). 

Therefore, Equation (35) can be reduced to 


% <'> *2 


S Q <“„> 

tt 


2 3 

M Cu> 
n n 


2 ”2 

Instead of q^ (t), Jones determined q^ , the mean square 

acceleration. Note that the Fourier transform of q is 

^n 


q n = (i*)\ 


Therefore , 


S q ( " ) ' ^ 


(ico) S- (w) 

'•n 


M H ( to) H (in) 
n n n 


(37) 


(38) 


(39) 


The result for q^ is (Ref. 21) 


1 n ^ \ 

8 M 2 ? % * 

n 


Let 


s r f 

<n 


(40) 


(41) 


where is Che freest ream dynamic pressure and is a 
nondlmensionai aerodynamic excitation parameter. It foLlows that 
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2/2 


-<J) 

c n 


1 / 2 ] 




M 


(42) 


(43) 


2/2 ' n 

To use Equation (43), the damping ratio ( £) is needed. It 

consists of both the structural damping (? s ) and the aerodynamic 
damping (C a )» The latter arises from the effective angle of attack 
due to wing vibration and is given by 


2M C oj q = 2q SK TT 11 
n a n n V 


(44) 


where K, the aerodynamic damping parameter, is a nondimensional 
parameter depending on the mode shape, the wing planform, and the 
sectional lift— curve slope. Equation (44) is assumed applicable to 
both attached and separated flows. It follows from Equation (44) 
that 


_ q ~ SK 

H u) V 
n n 


K 


M u) VC 
n n a 


q 


S 

CD 


(45) 


(46) 


In Jones method , both E and EC are assumed to be independent of 
the scale effect. In other words, their values determined from 
model test can be applied to full-scale airplanes. Practical 
procedures of applying this method were discussed by Butler and 
Spavins in Reference 24. They are as follows: 

(1) Determine modal frequency the mode shape, generalized 

mass Mjj , and structural damping from wind— off resonance 
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tests on model and aircraft. Note that the relevant model 
mode shape must be approximately correct. 

(2) Measure rms acceleration or bending moment C„ at a point 

D 

on the wing, the total damping flow velocity V, and 
dynamic pressure at a given Mach number and angle of 

attack in wind-tunnel tests. 

(3) Relate C to q in generalized coordinates using the mode 

B n 

shape (see Section 2.3). 

(4) Calculate E from Equation (42). 

(5) Calculate K from Equation (46). 

(6) Calculate total damping of aircraft by adding calculated 

5 from Equation (45) to the measured C . 
a s 

(7) Predict rms acceleration or bending moment at a point on 
the aircraft wing from Equation (43) using the measured 
aircraft mode shape. 

Mabey f s Method (Refs. 25 and 26): 

This method was developed to determine qualitatively the flight 
conditions for light, moderate and heavy buffeting for the full- 
scale aircraft from measurement of wing root bending moment of a 
conventional wind-tunnel model. It is assumed that the wing 
responds to buffeting pressures in somewhat the same way as to the 
wind-tunnel turbulence at the wing fundamental frequency. 

Let the tunnel unsteadiness /nF(n) be defined so that the total 


rms pressure fluctuation coefficient is given by 



where 


a = f W/V 
o 

W = tunnel width 

f = wing fundamental bending frequency in cycles 
per second 

V = freestreara velocity. 

Define 

C bb (M, a) = wing-root strain signal/q,* (48) 

Before the onset of flow separation on the model, C dd (M, a) has been 

BB 

shown experimentally to be constant equal to C (M, a = 0) . This is 

BB 

the portion of the model response caused by the tunnel unsteadiness 
/nF(n). Assume that 

c BB^ M » a = 0°) = Kg/nF(n) (49) 

where Kg is a scaling factor. Then 

C BB a * 0°) = a = 0°) = ^nF(n) 

B 

(50) 

Beyond buffet onset, Cgg(M,a) is increased due to wing buffet 
pressures. Let 

C bb "(M, a) = [C B8 '(M, a) 2 - Cg B '(M, a= 0°) 2 ] l/2 (51) 

The angLe of attack at which Cg B "(M, a) first differs from zero is 
the buffet-onset angle. From correlations on nine models of fighter 
aircraft, the following buffeting criteria were suggested: 



Buffet onset 


0 


Light 

buffeting 

c bb" 

= 0.004 

Moderate buffeting 

c bb" 

=• 0.008 

Heavy 

buffeting 

C BB 

= 0.016 


Note that in using this method, the total damping of the wing 
fundamental mode should be relatively constant, independent of wind 
velocity and density* This is true if models with solid wings of 
steel or light alloy are used, because in this case the structural 
damping will predominate. No mass, stiffness (or <a> 0 ) and damping 
for both models and aicraft are needed. It is useful during 
comparative tests for projects with alternative wing designs. 

Thomas' Method (Ref. 27): 

At transonic speeds, buffeting is closely connected with flow 
separation due to shock-boundary layer interaction and shock 
oscil Lat ions . Using conventional boundary layer methods, the 
development of boundary layer on airfoils at transonic speeds can be 
calculated. By comparing calculations with experimental results, 
Thomas postulated that buffet onset started if the point of rear 
separation coming from the trailing edge reached 90% of the airfoil 
chord. 

Redeker (Ref. 28) extended this method to infinite yawed wings 
by using the pressure distribution on a section normal to the 
leading edge and applying a three-dimensional compressible boundary 
layer method. 


Lb 



Further extension of Thomas’ method to finite wings was made by 
Proksch (Ref. 28). A buffeting coefficient (C Bi ) is defined which 
is directly related to the rms value of the wing root bending 
moment. It is assumed that the fluctuations of the wing root 
bending moment are proportional to the integral evaluated along the 
wing span of the product of local lift fluctuations and the distance 
from the wing root ( n - n^) • A further assumption is that the local 
lift oscillations caused by flow separation are proportional to 
length £ g (n) of the separated flow at a spanwise station of the 
wing. It follows that 

1 V n) ZZT 

c Bi = / — i — (n - n R )dn ~ /c (52) 

n R c 

2.3 The Present Proposed Method 

Theoretically, it is possible to use Equations (23) - (28) to 
calculate buffet response in the most general way. However, it 
would be an expensive undertaking because extensive fluctuating 
pressure measurement on empennage must be made. In addition, these 
fluctuating pressures are configuration dependent and vary with 
flight conditions. Therefore, a method similar to Jones' in 
structural representation is proposed. That is, the structural 
motion is assumed to be governed by a single— degree - of "freedom 
system for each mode (see Equation 32). However, Jones' method is 
used primarily to extrapolate model test results to full-scale 
aircraft. On the other hand, in the present method the buffeting 
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E 

excitation, and hence the generalized aerodynamic force Q (t), is 

n 

calculated from given configuration geometry, based on a limited 
amount of force measurements on delta wings. In this report, only 
the calculation of buffeting force will be considered. The final 
solution of Equation (32) and its applications In calculation of 
root bending moment will be Investigated in the future. It should 
be noted that the proposed method accounts only for buffeting due to 
vortex breakdown. 

In developing the proposed method in calculating the buffeting 
force, the following steps are needed. 

(1) Buffeting vortex strength in the burst region must be 
known. It is known that steady vortex strength from a 
slender wing or LEX can be estimated by the method of 
suction analogy (Ref. 29). Similarly, buffeting vortex 
strength can also be estimated If buffeting normal force 
data on slender wings are available. This is because any 
buffeting on slender wings can be assumed to be caused by 
the leading edge vortex. A limited amount of such data 
was published in References 4 and 9. 

Let c g be the sectional suction coefficient. Based 
on the suction analogy, the vortex lift is proportional to 
c g . The vortex lift can also be expressed in terms of the 
vortex strength V through Kutta-Joukowski theorem as 

I PV 0o “ c s cc, y = prw e e d£ ( 53) 
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where w^ 0 is the normal velocity at the leading edge and 
d£ is the vortex length along the leading edge. It 
follows that 


c c 
s 


V dZ 2 w /V 
00 £e 


and 


r t - /£-<** - t/ 


dy 

b/2 


c c 
s 


w- /V 
£e » 


dy 


The average strength per unit length is 


V t ' T t'\e 


- I 

Ze o 


b/2 


c s c 

W /V 
£e 00 


dy 


(54) 


where is the length of the leading edge. The unsteady 

aerodynamics program of Reference 30 was revised to 
calculate r * 

In addition, a line unsteady vortex from the LEX-wing junction 
is used to generate buffeting flow for vertical tails. This line 
vortex will produce additional loading, to be called "augmented 
vortex lift," due to momentum transfer. According to Figure 1(a), 
if the force due to momentum transfer is equated to the vortex lift 
through the suction analogy, then 


/ J P v £ c c dy = / P r w dl 


= / pV'(V' • dA) 

= / 2p(i V ) (1 V ) d£ 7 

v 2. 00 J 2 00 av 


(55) 
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where V f is taken to be -L V based on available data (Ref. 31). 

I °° 

From Equation (55), it follows that 

7 = I / C c dy (56) 

av S „ ; s 

£e 

The augmented vortex lift on a downstream lifting surface is 
then given by (see Figure 1(b)) 


Augmented Vortex Lift = 2 p(y V^J c r 


■ 1 ° V i 2 r t (f)*e ' (57 > 


In the calculation, the buffeting normal force Is obtained by 
assuming a vertical oscillation of constant amplitude over the 
region of predicted vortex breakdown. The latter was calculated by 
a serai-empirical method to interpolate or extrapolate experimental 
data (Ref. 32). The amplitude was adjusted to match the 
experimental data on mean square values of fluctuating normal force 
coefficients given in Reference 9. The resulting unsteady leading- 
edge suction is then used to calculate the buffeting vortex 
strength. Unfortunately, only data at a low frequency for some 
delta wings were measured in Reference 9. On the other hand, the 
power spectrum over a range of frequencies at the vortex-breakdown 
angle of attack for the BAG 22 1 configuration is available (Fig. 24 
of Ref. 4). This is illustrated in Figure 2. Unless additional 
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data are available in the future, for the present purpose the low- 
frequency data of Reference 9 will be used to derive the buffeting 
vortex strength for a range of angles of attack. At other 
frequencies, the power spectrum is assumed to be that for the BAC- 
221, and the strength will be multiplied by a ratio obtained from 
data for the BAC 221 In Reference 4, 

(2) RMS root bending moment can be calculated as 
b/2 n 

M o (t) = ^ + V y> “ ( 1 9 n ( t ) < t’ n (y)]ra(y)]ydy 

° n=l 

(58) 

where £ g is the sectional lift due to external forces, 
is the sectional lift due to structural motion and the 
last term is the inertial forces. For a rigid wing, the 
last two terms can be ignored. In Jones' analysis, was 
also ignored. Using the notation of Equation (4), fc can 
be written as 

N N q . 

£^( y > t) = q^ £ r~ / AC (x, y) <£ (x, y)mdx (59) 

n=l j=l o p j n 

Let the Fourier transform of £ g (y, t) be written as 

N 

l E (y ’ c) = I * E (y, iw)Q E (iu) (60) 

n=l j q 

where £ £ is the sectional lift due to a unit generalized 
n 

force in the n 1 -* 1 mode. Applying the Fourier transform to 
Equation (58), it is obtained that 
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V N 2 
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5 0 (i«) - / i n E <y, 1.3)+^+ i. «fo:)"+ 

o n 3 l n Q n^l a n 


I i bm e (1 “> + BM « < lB > + m-tuo /*+„ ( y ) y d y' q n 

n=l n a n n 


^ E 




( 61 ) 


where Equation (33) has been used- ^ 3 ^( a} ) the bending 
moment transfer function and is defined as 


N 2 b/2 

h BM ( “) ■ l l BM E < l “> + BM H <1u) + M-rM I "*„<>’»' dyI 
n=*l n n “ - 


n n o 


(62) 


The power spectral density of M (t) is therefore 


S BM (W) " I H BM (W) I S Q 

n 


(63) 


where S- is the power spectral density of the buffeting 
excitation. For a rigid wing, Equation (63) can be 
simplified to 


S BM ( C ° ) = i BM E (itAj) i S 0 ((AJ) 


(64) 


In applications, will be calculated by assuming a 

unit buffeting excitation over the region of vortex 
breakdown at a range of frequencies. The mean square 
value of root bending moment is then given by 


M o = / S BM (uj)dU) = 2 / S BM (u,)da) 

-oo O 


(65) 
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which is to be integrated numerically . Square root of the 


integrated value provides the rms value of root bending 
moment . 

(3) Since only total force power spectrum, instead of pressure 
power spectrum, will be used, it is assumed that the 
pressure fluctuations at every point on the wing are 
perfectly correlated in space and are in phase. Based on 
this assumption, Mabey and Butler showed that the total 
force power spectral density was proportional to the 
pressure power spectral density (Ref. 11). The results 
from this were shown to be reasonably accurate. 

In the present application to empennage buffeting due 
to a LEX vortex, those unstedy buffeting vortices, once 
generated, will be convected downstream in accordance with 
the general principle of unsteady aerodynamics. 

(4) With the power spectral density of buffeting vortex 
strength determined at a given flight condition, 
fluctuating normal velocity will be induced on the 
empennage. By satisfying the usual flow tangency 
condition, buffeting pressure spectral density on the 
empennage can be calculated. From the buffeting pressure 
spectral density, the power spectrum of bending moment or 
other aerodynamic characteristics can be determined. The 
mean square values of root bending moment are calculated 
by using Equation (6 5). 
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(5) Similar to Jones 1 method, the calculation of buffet 
response requires structural data, such as mode shapes, 
generalized mass, and damping ratio- Aerodynamically 
induced damping can have major effect on buffet response 
and it must be accounted for in analytical representation 
of flexible aircraft. 

(6) In applications to empennage buffeting, the locations of 
LEX vortex bursting may be based on experimental data or 
theoretical calculation (Ref. 32). 
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3. 


NUMERICAL RESULTS 


Extensive measurements of buffeting normal force were made by 

/ 

Earnshaw and Lawford on a series of delta wings at a low frequency 
of n = 0.05, where 


n = fc/V 

oo 


( 66 ) 


and f is the frequency in cycles per second. "n“ is converted into 

the conventional reduced frequency k by multiplying by 2 tt. Before 

vortex breakdown, the normal force fluctuation is assumed to be 

caused by tunnel flow unsteadiness. At each angle of attack, a 

buffeting vortex strength T^ can be calculated from Equation (54). 

The same expression is used to calculate the vortex strength "f in 

s 

steady flow using the steady-flow c g . if the results are plotted in 
terms of the ratio of buffeting to steady vortex strengths, , : 

R b/s = V T s (67) 

curves for cambered and flat wings tend to collapse into a single 
one. This is illustrated in Figure 3 for a 70-degree delta wing. 

In Figure 3, Aa is the incremental angle of attack beyond that of 
vortex breakdown at the trailing edge. It follows that the 

buffeting vortex strength is a function of steady-flow vortex 
strength and Aa. 


To extend limited available data to different planforms, a 
correlation parameter capable of correlating vortex-breakdown 
characteristics is needed. In Reference 32, the nondimensional 
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distance measured along the leading edge from apex to the centroid 
of leading-edge suction distribution, y^ , was found to be a useful 
parameter for this purpose. If this y^ is used to analyze all 
experimental data for a series of delta wings in Reference 9, 
results for buffeting vortex strength can be plotted in one graph as 
shown in Figure 4. Note that the dash curves in Figure 4 represent 
extrapolation, because data in that range of y ^-values (i.e. highly 
swept delta wings or LEX) are not available. 

3.1 Results for a 65-degree Delta Wing 

To check the theory, test data of Reference 33 for the root 
bending moment of a 65-degree delta wing will be used. Static 
bending moment coefficients based on c are presented in Figure 5. 
Calculated results from Reference 29 are also presented for 
comparison. It is seen that at high angles of attack, the theory 
overpredicts the root bending moment. The predicted is about 20 

degrees, while it is about 18 degrees in the experiment. 

To find the buffeting characteristics, it is assumed that the 
buffeting excitation (i.e. the plunging amplitude in the vortex- 
breakdown region) for the flat wing is the same as that for a 
cambered 65-degree delta wing if R^/ s Is the same. Therefore at a 
given Aa, is obtained from Figure 6. Using this R^/g, the 

buffeting excitation (AMPLG) can be determined from Figure 7. Note 
that Figure 7 was constructed from the experimental data of 


26 



Reference 9 for a cambered wing. The resulting buffeting excitation 
amplitudes for a flat 65-degree delta wing are plotted in Figure 8. 

To calculate the dynamic response of a rigid wing. Equations 
(64) and (65) will be used. To calculate the transfer function 
(BMg) for the root bending moment, a unit amplitude of vertical 
excitation is prescribed over the region of vortex breakdown at each 
frequency. Some results are presented in Figure 9 • The 
corresponding power spectral densities for the excitation are 
obtained by multiplying the values in Figure 8 (for a low frequency 
only) by a ratio obtained from Figure 2 for other frequencies. The 
results are shown in Figure 10. Equation (65) is then integrated by 
the trapezoidal rule to produce the mean square values of root 
bending moment. The rms values are presented in Figure 11. Note 
that experimental data were obtained at resonant frequencies of the 
fundamental bending mode. Since the spectral density is higher at 
higher frequencies (Fig. 10), the calculated response of a rigid 
wing tends to be similar to the test data at a high frequency, 
although the magnitudes are underpred icted . It is expected that the 
prediction can be improved if the structural flexibility is 
accounted for. 

3.2 Results for an F-18 Configuration 

Aerodynamics calculation for an F-18 conf iguration is based on 
the code of Reference 29* The modeling of geometry is shown in 
Figure 12. In the calculation, wing sectional aerodynamic 
characteristics are needed to account for the effect of viscous 



separation at high angles of attack. These characteristics are 
calculated with Eppler's code (Ref. 34). These calculated results 
are modified for a beyond the stall a by using experimental data for 
a flat plate. The strake (i.e. LEX) vortex position and its 
bursting point can also be calculated and are correlated with water 
tunnel results in Figure 13. In this calculation, the longitudinal 
location of vortex bursting point at a = 30 degrees is adjusted to 
match that from water tunnel results. At other angles of attack, it 
is determined by extrapolation based on data on bursting progression 
rates for delta wings (Ref. 32). Additional water tunnel results 
are presented in the Appendix. The predicted longitudinal 
characteristics are presented in Figure 14. It is seen that the 
results are reasonably well predicted. 

The calculated LEX vortex position and bursting points at 
different a’s are then used in the unsteady aerodynamics program to 
calculate the fin buffeting. The buffeting vortex strength is 
determined from Figure 4, In the calculation, effects of both 
induced normal velocity from the buffeting vortex and augmented 
vortex lift are included. The resulting rms values of root bending 
moment without structural flexibility are presented in Figure 15. 
Since there are no appropriate test data for comparison, application 
of Figure 14 is illustrated with the following numerical example. 
Assume that ot = 2 5 degrees and = 30 psf. From Figure 15, C 

Bd 

is found to be 0.00765. Since the reference length is 11.12 ft. 

( =b/2) and the reference area is 104 ft^ (= S re f) for the vertical 
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tails, it follows that the root mean square value of root bending 
moment (= BM) is 


BM 


^Bd ^ref 


b/2 


= 0.00765 x 30 x 104 x 11.12 = 265.4 ft 


lb 
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4. CONCLUSIONS AND RECOMMENDATIONS 


Test data showed that vortex breakdown could cause severe 
buffeting on a lifting surface- In the present Investigation, a 
method was developed to describe buffeting flow after vortex 
breakdown by interpolating available data through a correlation 
parameter. The buffeting flow is used as the aerodynamic forcing 
function to calculate surface pressure distribution and root bending 
moment through an unsteady lifting-surface theory. 

To extend the present work, the following are recommended. 

1. The buffeting flow field of a vortex be measured at some axial 
and circumferential locations. A theoretical method based on 
statistical concepts, such as those used In inlet dynamic 
distortion (Refs. 36 and 37), should be developed to describe 
the buffeting flow field. The latter can be assumed to be 
Independent of tail geometry and is used to generate tail 
pressure distribution through a lifting surface theory. This 
approach Is different from the conventional method of buffet 
prediction which requires the measurement of tail surface 
pressure distribution directly. 

2. The simplified aeroelastic model equation (Eq. 32) be solved 
with the forcing function provided under Task 1. The bending 
moment is calculated through Equation 58. The combined results 
of Tasks 1 and 2 should provide a quick estimate of buffeting 
level on vertical tails when geometry of the latter is changed. 
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Buffeting Excitation for a 65-deg Flat 
Delta Wing. n=0.05 % 
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Figure 9. Root Bending Moment Transfer 
Function for a 65-deg delta 
Wing . 
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Figure 10. Power Spectral Density of 
Buffeting Excitation for a 
65-deg delta Wing. 
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Figure 12. F-18 Geometry Modeling and Strake Vortex Location 

at alpha=30-deg . 
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WATER TUNNEL DATA 
VORSTAB CALCULATION (Ref. 29) 
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Breakdown for an F-18 Configuration 


WATER TUNNEL DATA 
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VPI DATA M-0.0 (Ref. 35) 
NAVY SIMULATION DATA M-0.2 
VORSTAB CALCULATION M->0.2 
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Figure 14 . Longitudinal Aerodynamic Characteristics 
for an F— 18 Configuration 


48 





Figure 15. Calculated RMS Root Bending Moment Coefficient for the 

F-18 Vertical Tail at M-0.0, Rigid Configuration, C Bd is 
Based on b/2 of the Vertical Tail Rotated to the X-Y Plane 
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ABSTRACT 


Water tunnel tests have been conducted to study the 
flow associated with fin buffet for twin-fin fighter air- 
craft using l/48th scale F/A-18 models. Flow visualization 
made use of colored dyes to determine vortex patterns, and 
surface hot film anemometry was used to study the turbulent 
energy and the frequencies present in the flow. Configur- 
ations tested included the full airplane, airplane without 
fins, airplane without leading-edge-extensions, (LEX's) and 
airplane without wings. Test Reynolds number ranged from 
4,300 to 12,800, with corresponding Mach numbers less than 
10 -6 . 

The flow studies show that the LEX vortices burst just 
forward of the fins at about 25° angle of attack. Removing 
the fins had negligible effect on vortex locations and 
bursting, but removing the wing had a marked effect on both 
location and burst angle of attack for the vortices. Stud- 
ies of body vortices with the LEX's removed demonstrated 
that the body vortices were not a dominant feature of the 
flow associated with fin buffet. 

Hot film anemometer signals show that fin surface tur- 
bulence increases with angle of attack, and that dominant 
frequencies appear in the flow when bursting occurs. The 
dominant frequencies correspond to a Strouhal number of 0.7 
for all speeds tested, and for all angles of attack for 
which vortex bursting was present. Flow patterns, vortex 
bursting angles of attack, and Strouhal numbers of the 
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unsteady flow correlate well with wind tunnel tests from 
other investigations at higher Reynolds number, confirming 
the validity of water tunnel testing. 

SYMBOLS 

c wing mean aerodynamic chord, (0.24 ft, l/48th 
model scale) 

f frequency, Hz 

M Mach number, V/ (speed of sound), non-dimensional 

Psd power spectral density of hot film signal, 

(Volts) 2 /Hz 

Re Reynolds number, cV/v , non-dimensional 

S Strouhal number, fc/V, non-dimensional 

V tunnel velocity, ft/sec 

a angle of attack, degrees 

v kinematic viscosity, ft^/sec 

INTRODUCTION 

Twin-fin arrangements have recently emerged as a con- 
figration favored by aircraft designers. This configuration 
is especially attractive for carrier-based aircraft, since 
it offers reduced fin height, making hanger access and 
maintenance easier. The F/A-18 aircraft uses this arrange- 
ment, but unfortunately the aircraft has developed fin 
fatigue problems requiring structural modification. Flight 
and wind tunnel tests revealed that the un-anticipated fin 
loads occur at subsonic high angle of attack conditions. 
These loads are apparently related to the interaction of 
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vortices emanating from the wing leading edge extensions 
(LEX's) with the fins. 

This report documents tests using small scale models in 
a water tunnel to visualize the flow phenomena associated 
with high angle of attack conditions for the twin-fin 
fighter type aircraft. The primary purpose of this research 
was to identify the flow associated with fin buffet for this 
aircraft and to generalize so far as possible from these 
results, in order to avoid such buffet problems for future 
designs. A second purpose was to evaluate the water tunnel 
as a research tool as compared to the more traditional wind 
tunnel and flight test environments for experimental tests. 

EXPERIMENTAL FACILITY 

The facility for the experiments reported here was the 
NASA Ames-Dryden water tunnel flow visualization facility 
located at Edwards Air Force Base, California. This tunnel 
is a closed return vertical flow water tunnel, with 16" x 
24" test section. Earlier tests of the F/A-18 in a similar 
water tunnel were reported in reference 1. The tunnel was 
designed primarily for use as a visualization facility, but 
in the present tests special surface hot-film anemometry 
instrumentation was utilized to make quantitative measure- 
ments of the unsteady, buffeting flow. 

LIMITATIONS OF SMALL-SCALE AERODYNAMIC TESTING 

Classical design of fluid dynamic experiments requires 
"dynamic similarity" of the model and full-scale airplane. 
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Dynamic similarity is achieved when Reynolds number and Mach 
number of the model and full-scale airplane are the same, 
and when the model and full-scale are geometrically similar. 
In practice, geometric similarity is nearly always achieved 
by using properly proportioned models. Experiments and 
theory have shown that matching Mach number is necessary 
only when compressibility effects become important. This is 
typically at Mach number above 0.6, depending on thickness 
ratio. For higher Mach numbers, the pressure distributions 
are directly affected by compressibility, and Mach matching 
is essential. 

For measurement of skin friction, and precise matching 
of separation and stalling of airfoils, matching Reynolds 
number is required. While Reynolds number matching is 
required in principle, in practice small-scale testing is 
frequently used; even though it almost always results in 
Reynolds number below full-scale values. Full-scale, 
pressurized and cryogenic wind tunnels are facilities in 
which full-scale Reynolds number is ordinarily achieved. 
Testing in these facilities is very expensive because of 
model and operational costs. Fortunately, it is the nature 
of viscous flow that aerodynamic characteristics are rela- 
tively insensitive to Reynolds number. Often Reynolds 
number differences of factors of 3 or even 10 have relative- 
ly small effect on all aerodynamic coefficients except 
parasite drag and maximum lift coefficients. For the 
particular case of fighter-type aircraft, which are charact- 
erized by thin, highly swept surfaces, operation at low 
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speed invariably implies high angles of attack, and high 
angle of attack leads to separation along the leading edges. 
For such cases, the separation locations are fixed by 
geometry, and aerodynamic force coefficients and pressure 
coefficients are essentially independent of Reynolds number. 
This is borne out by the research of reference 2. It is 
this peculiar combination of sharp leading edges and high 
angles of attack that lends validity to small-scale smoke 
tunnel and water tunnel tests of aircraft designed for 
supersonic flight. Test planning and interpretation of 
results must recognize the limitations and the regimes for 
which testing at Reynolds number substantially lower than 
full scale may yield useful information. 

MODELS 

Models were variations of the F/A-18 aircraft, fabri- 
cated from 1/48 scale hobby shop kits. The hobby kits are 
dimensionally sufficiently accurate for tests of this type, 
so the only modifications required were the addition small- 
bore tubing to accommodate the dye used for stream tracing, 
and a mounting strut. Dye tubes were connected to manifolds 
within the models, which were fed from a separate external 
dye reservoir for each of the colors desired. In addition, 
the models were equipped for engine inlet simulation by use 
of hollow nacelle passages and flexible plastic water siphon 
tubes attached to the engine exhaust nozzles. By drawing 
water into the inlets through the siphon tubes at an appro- 
priate rate, it is possible to control the engine mass flow 
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capture area ratio. This permits proper simulation of 
streamline patterns around the aircraft. This is partic- 
ularly important for those regions of flow near the engine 
inlets, such as the wing leading edges. 

Since one purpose of these tests was to evaluate the 
effects of the various aircraft components on fin buffet, 
models were designed with several geometric variables. It 
was less expensive to construct a family of four models with 
the various configurations than to construct removable hard- 
ware for a single model. Wings leading flaps were deflected 
34 degrees down, and trailing edge flaps were undeflected 
for all tests. These settings are consistent with flight at 
angles of attack 25 degrees and greater. The geometric 
variables tested are described below. 

F/A-18 BASIC MODEL - Complete airplane with leading 
edge flaps deflected 34 degrees, trailing edge flaps 
neutral, all tail surfaces neutral. 

F/A-18 WITHOUT WINGS - Same as basic model, except 
wings removed outboard of leading edge extensions (LEX's). 
This model was used to evaluate the interference effects of 
the wing and flow fields. 

F/A— 18 WITHOUT FINS — This model was used to evaluate 
the possibility that the fin "blockage" might generate an 
adverse pressure field of sufficient strength to cause pre- 
mature bursting of the leading edge vortices. 

F/A-18 WITHOUT LEX'S - The purpose of this model was to 
identify the role and interaction of forebody and LEX vor- 
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tices, and to ascertain possible wing or forebody vortex 
interactions on the fin. 

TEST CONDITIONS 

For all tests, inlet flow was established to provide 
for an inlet capture area ratio of unity. Speed control in 
the water tunnel is by means of a valve with a series of 
fixed settings, preventing infinitely variable speed con- 
trol. Table 1 gives the speeds used for these tests, and 
corresponding chord Reynolds numbers and Mach numbers. 


Table 1 - Test Conditions 



Reynolds Number 

Mach Number 

0.25 

4,300 

0.5 x 10~ 6 

0.58 

9,900 

1.3 x 10 -6 

0.75 

12,800 

1.6 x 10 -6 


Angle of attack was varied from 0 to 40 degrees, in 
increments of 5 degrees. At 40 degrees, the model was 
nearly in contact with the upper wall, so higher angles 
could not have been tested without the use of an offset 
sting mount. Video and still pictures were obtained from 
top and side views in separate runs. 
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INSTRUMENTATION 


Instrumentation consisted of video cassette recording 
equipment and conventional camera for still photos. In 
addition, the fin of the basic model was fitted with a 
strain-gage at the fin root, and two different types of 
surface hot film anemometers (Disa and Micro-Measurements). 
These instruments were intended to detect flow unsteadiness 
over the fin, for correlation with dynamic strain gage data 
from flight tests. 

The Disa hot-film anemometer and the strain gage pro- 
vided very low-level signals, and these could not be distin- 
guished from random background noise. Only the Micro- 
Measurements hot film gage provided a signal which displayed 
characteristics which changed in a consistent manner with 
angle of attack. Therefore, only the data from the Micro- 
Measurements hot-film gage was utilized for dynamic measure- 
ments. This gage was located at 63% span and 50% chord on 
the inboard surface of the starboard fin. The output signal 
from this sensor was monitored on an oscilloscope, and 
selected signals were also processed using a modal analyzer. 
The oscilloscope provided real time characteristics of the 
signal, and the modal analyzer provided frequency analysis 
of the data in the form of power spectral density (Psd) 
graphs of the gage voltage. It should be noted that the hot 
film gage utilized in this manner provides a measure of the 
heat transfer at the surface. This signal provides a quali- 
tative but not quantitative measure of surface skin 
friction. The merit in this instrument is the ability to 
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extract information about the frequencies present in the 
unsteady flow. Futhermore, these signals can be processed 
to obtain qualitative differences in turbulent energy levels 
for the various angles of attack and tunnel speeds. 

RESULTS OF FLOW VISUALIZATION 
BASIC F/A-18 - The flow video and still pictures show a 
consistent and repeatable pattern for the vortex flow of 
this aircraft configuration. Vortex "bursting" or "break- 
down" is characterized by an abrupt increase in vortex core 
diameter, often preceded by a spiralling of the core before 
complete turbulence ensues. All vortex burst patterns 
exhibit some unsteadiness, with burst locations oscillating 
somewhat with time. For this reason, the video images were 
used as the primary source to determine "average" burst 
locations. Geometry and vortex burst locations are shown in 
figure 1. As angle of attack is increased from 0°, vortices 
form along the LEX's. These vortices increase in strength 
with angle of attack, and flow aft above the horizontal tail 
surfaces but beneath the fins for angles below 20°. At 20°, 
vortex bursting occurs aft of the wing trailing edge and 
outboard and beneath the fin. At 25° angle of attack, the 
LEX vortex burst point is located more forward and inboard, 
with the axis of rotation nearly in line with the fin 
leading edge. The burst point is slightly forward from the 
fin leading-edge at this angle of attack. In figure 1 (b), 
sketches of vortex flow field after bursting are presented 
to show the proximity of this flow field to the fins. It is 
difficult to conceive of a vortex burst location which could 
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be potentially more detrimental in terms of introducing fatigue- 

producing loads. As angle of attack is increased further, 
the burst point moves progressively further forward. 

F/A-18 WITHOUT FINS - The absence of the fins had 
negligible effect on the flow field. Vortex locations and 
burst positions were indistinguishable from the basic model. 
These results indicate that the pressure gradients associ- 
ated with the fins are not strong enough to instigate vortex 
bursting. This is not surprising, considering the small 
thickness-to-chord ratio of the fins. 

F/A-18 WITHOUT WINGS - In this configuration, initial 
formation of the LEX vortices was in much the same manner as 
for the full model. As angle of attack was increased, 
however, the vortices were located more inboard than on the 
basic model, and they remained intact, without bursting, up 
to 30° angle of attack. This test series shows important 
changes in vortex locations and a substantial delay in 
vortex bursting when the wing is removed. When the wing is 
present, the adverse pressure gradient field over the aft 
portion of the wing is evidently a dominant factor in pro- 
ducing vortex bursting. 

F/A-18 WITHOUT LEX'S - At high angles of attack, the 
fuselage forebody, like the LEX's, produces a pair of vor- 
tices. This test series shows that the fuselage vortices 
initially form inboard of the fins, and that they remain 
inboard of the fins over the entire angle of attack range 
for which they are visible. At 15° these vortices trail 
between the fins without bursting. At 20°, bursting occurs 
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slightly aft of the wing trailing edge, and at 25° the burst 
location is near wing mid-chord. Burst location remains 
essentially unchanged as alpha is increased from 25° to 40°. 
Even at 40°, the spiraling wake of the burst vortices 
remains close to the fuselage centerline, inboard of the 
fins. This test series shows that the forebody vortices 
have much smaller interaction with the fins than the LEX 
vortices. Studies of the video tapes also show that the 
spiraling of forebody vortices is at a distinctly lower rate 
than LEX vortices, indicating lower vorticity associated 
with the forebody vortices, in agreement with theory. 

SUMMARY OF FLOW VISUALIZATION TESTS - Results of the 
flow visualization test series show that the LEX vortices 
are the dominant flow feature which provides strong inter- 
action with the fins, and that this interaction is maximum 
in the 25° to 30° angle of attack range. In this angle of 
attack range, bursting occurs just ahead of the fin leading 
edge. Bursting is not influenced by the fins, but is 
strongly influenced by the presence of the wing. 

RESULTS OF HOT-FILM SIGNALS (BASIC MODEL ONLY) 

Recently, the "modal analyzer" has been developed for 
interpretation of dynamic test data, particularly from 
vibration and flutter testing. This device provides high 
rate analog-to-digital conversion of signals, with digital 
storage and processing, including fast Fourier transform 
techniques (FFT) for determining frequency content of a 
signal. In addition to its application to strain gage and 
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accelerometer signals, the modal analyzer is well-suited for 
analysis of hot film anemometry signals from unsteady or 
turbulent flows. In the case of fin buffet induced by 
vortex interactions, it was anticipated that the vortex 
impingement or vortex bursting process might be periodic, or 
at least have a characteristic frequency signature which 
could be used for vibration and fatigue analysis. 

The modal analyzer was used to obtain power spectral 
density (Psd) data for each angle of attack from 0° to 40°. 
Integration of the Psd data over the frequency spectrum 
yields energy content of the fluctuating voltage across the 
hot film gage. The Psd integral is, in turn, a measure of 
the turbulent energy content of the airstream near the gage. 
Figure 2 is a graph of energy versus angle of attack. While 
these data show considerable scatter, the figure clearly 
indicates a trend of increasing energy content with angle of 
attack. 

Figures 3 through 7 show the power spectra for the 
various angles of attack and speeds. Peaks in the Psd plots 
indicate frequencies which are characteristic of the flow. 
Thus a Psd graph which is relatively flat indicates a flow 
with no identifiable periodicity ("broad-band" turbulence). 
In contrast, a Psd graph with distinct peaks is indication 
of characteristic periodicity within the flow. For angles 
of attack of 0° to 20°, Psd's are very low-level over all 
frequencies, and without peaks to indicate a dominant fre- 
quency. Since all Psd plots from 0° to 20° appear the same, 
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only the 20° graph is shown (figure 3). For 25°, a distinct 
peak appears at V = 0.25 ft/sec, but is not visible for V = 
0.58 and 0.75 ft/sec. For angles of attack of 30°, 35° and 
40°, energy levels are distinctly higher, and dominant 
frequencies are discernable at all speeds tested. 

To correlate model and full scale periodic phenomena or 
periodic phenomena for the same model at different speeds, a 
non-dimensional form of frequency is needed. This is pro- 
vided by the Strouhal number, which is defined as follows: 

Strouhal number = (f x c)/V, 

where: f = frequency, Hz 

c = wing reference chord 
V = free stream velocity. 

Presumably, Strouhal number, like lift coefficient and other 
aerodynamic coefficients, will be relatively invariant with 
speed and scale. To test for consistency of Strouhal 
number, runs were made with tunnel speeds of approximately 
two and three times the nominal value. Dominant frequencies 
selected from the Psd graphs are plotted as frequency versus 
velocity for each angle of attack from 20° to 40° in figures 
3 through 6. These results, taken from two separate test 
series, show that dominant frequency tends to increase 
linearly with tunnel velocity, indicating that Strouhal 
number is indeed constant with velocity. Furthermore, the 
Strouhal number is essentially a constant value of 0.7 for 
all angles of attack for which periodic behavior was 
observed. 
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CORRELATION WITH OTHER TESTS 


The observation that Strouhal number is approximately 
constant with velocity is a crucial finding from these 
experiments. These results lend validity to the use of 
small scale tests, since they show that Strouhal number is 
independent of Reynolds number, at least for the test range 
of velocities. Furthermore, wind tunnel tests of a 12% 
scale basic F/A-18 configuration at Reynolds numbers of 0.4 
and 0.8 x 10 6 by McDonnell-Aircraf t Company (ref. 3) have 
shown nearly the same location for LEX vortices, nearly the 
same angle of attack for vortex bursting, and nearly the 
same Strouhal number for unsteady fin surface pressures. 
These correlations are strong evidence that the fundamental 
flow patterns are independent of Reynolds number, even for 
the very low Reynolds numbers of the water tunnel tests. 
Achieving this correlation for the basic model gives cred- 
ibility to conclusions from water tunnel test results with 
the non-standard configurations. The merit of water tunnel 
experiments is the ability to quickly and inexpensively 
evaluate trends and essential features of flows associated 
with a wide variety of configurations. More detailed tests 
in wind tunnels and flight with a narrower range of config- 
urations are still essential to validate final design 
configurations. 
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CONCLUSIONS 


1. Surface hot-film anemometry shows high turbulent 
activity on the fins at conditions coincident with 
vortex bursting observed from flow visualization. 

2. LEX vortex bursting occurs directly forward of the fin 
leading-edge for angles of attack of 25° and higher. 

3. The onset of vortex bursting produces flow unsteadiness 
with a dominant frequency at a Strouhal number of 0.7 
for three speeds and for angles of attack from 25° to 
40°. 

4. LEX vortex bursting is associated with wing separation 
and stalling. Removing the wing produced substantial 
changes in vortex positions and delayed vortex bursting. 

5. Vortex bursting is unaffected by the fins. Removal of 
the fins had no appreciable effect on vortex locations 
or vortex bursting. 

6. Vortex frequencies, vortex bursting, and dominant fre- 
quencies from the water tunnel tests correlate well 
with wind tunnel tests at higher Reynolds numbers. 
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(a) Basic Model 

Figure 1. Configurations and Vortex Burst Locations 
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(b) Sketches of Vortex Flow Field After Bursting 
at an Angle of Attack of 25 Degrees 

Figure 1. Continued. 
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(c ) Model Without Wings 
Figure 1. Continued. 
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(d) Model Without LEX'S 
Figure 1. Concluded. 
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Figure 5. Power Spectral Densities, Alpha « 30° 


73 




volts /HZ volts /HZ voltsVHZ 










Frequency (Hz) Frequency (Hz) 


3.0 — 



(c) Alpha * 35° 



(d) Alpha - 40° 
Figure 8. Concluded. 


77 


1. Report No. 

NASA CR-1 79426 


2. Government Accession No. 


4. Title and Subtitla 

Investigation of Empennage Buffeting 


7. Author(s) 

C. Edward Lan and I.G. Lee 


9. Performing Organization Name and Address 


3. Recipient's Catalog No. 


5. Report Date 

March 1987 

6. Performing Organization Code 


8. Performing Organization Report No. 
H-1393 

10 Work Unit No. 

RTOP 505-68-71 


Flight Research Laboratory 

University of Kansas Center for Research, Inc. 

11. Contract or Grant No. 

Lawrence, Kansas 66045 

NAG2-371 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Contractor Report-Final 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

1 15. Supplementary Notes j 


NASA Technical Monitor: David F. Fisher, Ames Research Center, Dryden Flight Research Facility, 
Edvards, California 93523-5000, NADC Technical Monitor: Ed Desks, Code 6041, Naval Air Devel- 
opment Center, Warminister, PA 18974-5000. 

| 16. Abstract 


Theoretical methods of predicting aircraft buffeting are reviewed. For 
the buffeting due to leading-edge vortex breakdown, a method is developed 
to convert test data of mean square values of fluctuating normal force to 
buffeting vortex strength through an unsteady lifting-surface theory and 
unsteady suction analogy. The resulting buffeting vortex from the leading- 
edge extension of an F-18 configuration is used to generate a fluctuating 
flow field which produces unsteady pressure distribution on the vertical 
tails. The root mean square values of root bending moment on the vertical 
tails are calculated for a rigid configuration. 

Results from a flow visualization and hot film study in a water tunnel 
facility using a 1/48 scale model of an F-18 are included in an appendix by 
Dr. William H. Wentz, Jr. of Wichita State University. The results confirm 
that the LEX vortex is the dominant forcing function of fin buffet at high 
angles of attack. 


17. Key Words (Suggested by Author(s)) 


18. Distribution Statement 


Buffeting Unclassified — Unlimited 

Vortex breakdown 
Vortex flow 
Water tunnel 


Subject category 39 


19. Security Classif. (of this report! 

20. Security Classif. (of this pa gel 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

83 

A05 


•For sale by the National Technical Information Service/ Springfield, Virginia 22161 


























